ad30c817-4a14-4937-9388-f8bba8bbd5c0 3 3 12 decimal

&[DATE] &[TIME] - &[FILENAME]

t0 1 time :

Continued on length of support points 'j' j K X x Thiele Optimiz NONE j <= rows(k) N j 1 + el K j 1 + el : c j el N j el : a j el K j el : jj j 1 - (0 range c jj el a jj 1 + el : a jj el K jj el x X jj el - c jj el / + : 21 line for a 1 el 71 line :

Optimiz NONE ... Plug the RHS in Cf(x) below 6 K X x Thiele K 1 el x X 1 el - K 2 el x X 2 el - K 3 el x X 3 el - K 4 el x X 4 el - K 5 el x X 5 el - K 6 el / + / + / + / + / +

X Y Δ Thiele ε 1018 -^: n X rows : i 1 n range M i 1 el Y i el : for j 2 n range i j n range dM M i j 1 - el M j 1 - j 1 - el - : dM dM abs ε < ε M i j el X i el X j 1 - el - dM / : if : 21 line for for i 1 n range c i 1 el X i el : c i 2 el M i i el : 21 line for Δ 0 ≡ c M if 61 line :

XY 0 1.0475 0.25 0.4963 0.5 0.2755 0.75 0.1703 1 0.1136 1.5 0.0589 2 0.035 3 0.0157 4 0.0086 6 0.0035 102 mat :

X XY 1 col : Y XY 2 col : 21 line

X length 10

XK X Y 0 Thiele :

K XK 2 col : X XK 1 col : 21 line

x Cf K 1 el x X 1 el - K 2 el x X 2 el - K 3 el x X 3 el - K 4 el x X 4 el - K 5 el x X 5 el - K 6 el x X 6 el - K 7 el x X 7 el - K 8 el x X 8 el - K 9 el x X 9 el - K 10 el / + / + / + / + / + / + / + / + / + :

L Δ N inc X 0 : i 1 N range X i 1 + el L Δ i * + : X 1 el L : 21 line for X 31 line :

range/populate L X 1 el : Δ 0.001 : N 5000 : 13 mat U L Δ N inc : 21 sys

U U rows el 5

t0 1 time :

data i 1 U rows range vx i el U i el : vy i el U i el Cf eval : 21 line for vx vy augment 21 line :

data x Cf 2 1 sys

1 time t0 - 15 s *

t0 1 time :

data i 1 U rows range vx i el U i el : vy i el U i el Cf : 21 line for vx vy augment 21 line :

data x Cf 2 1 sys

1 time t0 - 42 s *